Entanglement Routers Using Macroscopic Singlets 
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We propose a mechanism where high entanglement between very distant boundary spins is generated by 
suddenly connecting two long Kondo spin chains. We show that this procedure provides an efficient way to route 
entanglement between multiple distant sites. We observe that the key features of the entanglement dynamics of 
the composite spin chain are remarkably well described using a simple model of two singlets, each formed by 
two spins. The proposed entanglement routing mechanism is a footprint of the emergence of a Kondo cloud in 
a Kondo system and can be engineered and observed in varied physical settings. 



Introduction:- A high entanglement between two well sep- 
arated qubits is the central resource for quantum communica- 
tion tasks. Entanglement between arbitrary pairs of distant 
qubits in a multi-qubit network enables the linking of sev- 
eral quantum registers to a single larger computer. It also 
facilitates the preparation multi-particle entangled states [l|], 
for measurement based quantum computation. One could 
ask whether many-body systems can serve as mediums for 
entanglement between arbitrary distant qubits in a multi-site 
network. Though this is the most important question from 
an "applied" perspective, the thriving field of entanglement 
in many-body systems |2J] remains focussed on the entangle- 
ment of blocks and proximal spins. In fact, long range entan- 
glement between individual spins is notoriously uncommon 
]3[]. There are proposals exploiting weak couplings of distant 
spins to a spin chain J4J, |5|], but these have limited thermal 
stability or a very long time-scale of entanglement genera- 
tion. Alternatively, a global quench [g] or specific time de- 
pendent couplings |7| may generate entanglement, though this 
decays with the system size. Finally, there is a proposal [8[] 
for distance independent entanglement through a local quench 
which, however, lacks the versatility of routing entanglement 
between multiple sites. A few quantum routers have been re- 
cently proposed 19|], but harnessing a canonical many -body 
phenomenon for routing still is an open question. 

Kondo systems IU0Ml2ll are very distinctive in the context 
of entanglement for at least two reasons. Despite being "gap- 
less", they support the emergence of a length scale £- the so 
called Kondo screening length IllOl lllll - which can be tuned 
by varying only one parameter [( 1 0f | and reflects in the entan- 
glement [8<], making it markedly different from other conven- 
tional gapless models. Furthermore, in Kondo systems, the 
impurity spin is maximally entangled 11211 with a block of 
spins whose spatial extent may be varied at will by tuning £. 

In this letter, we propose a dynamical mechanism by which 
long range distance independent entanglement may be gener- 
ated by the switch on of a single coupling suddenly connecting 
two macroscopic singlets. We show that this mechanism pro- 
vides an efficient way to route entanglement between various 
distant parties. By a macroscopic singlet we mean an arbitrar- 
ily long spin chain which has been engineered to behave as a 
Kondo system of pertinent £ and thereby as a two spin singlet. 



Indeed we show that, the key features of our mechanism for 
the two macroscopic singlets, are remarkably well described 
by a four spin system made of two singlets. 

Simple example:- Let us first consider two spin singlets 
each formed by only two spins interacting with a Heisen- 
berg interaction of strength J{ and J' 2 , respectively. The 
ground state of the composite system is then given by \gs) = 
\ip-) ® \ip~) with \ip~) = (|01) - |10))/\/2. In this sim- 
ple setting, one may generate high entanglement between the 
boundary spins, 1 and 4 by merely turning on an interac- 
tion J m between the spins 2 and 3. After quenching, the 
evolution of the system is ruled by the Hamiltonian H = 
J[ i?i.i72 + J2<73.<74 + J m o t 2-az and, since the initial state is 
a global singlet, time evolution allows for a nonzero overlap 
only with the singlet subspace of the spectrum of H so that 
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where, \S\) and IS2) are two singlet eigenvectors of H with 
energy Es 1 = — 4 J m and E$ 2 = respectively. In order to 
get maximal entanglement between the boundary spins 1 and 
4- after a certain time t*- one has to choose J m = J[ + ,/,. 
Once this condition is satisfied the state of the system at time 
t, up to a global phase, is given by 
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Surprisingly, \ip(t)} depends only on J m and, by tracing out 
the spins 2 and 3, one gets the density matrix pu(t) of the 
boundary spins. The entanglement between the spins 1 and 4 
may be easily computed using concurrence 11311 yielding 



l-3cos(4J m i) 
h = ?7iaa;{0, }. 



(3) 



Eq. d3j shows that E oscillates with a period of -^j— and that, 
at time t* = jj— , the spins 1 and 4 form a singlet state. In this 
simple setting one sees that: (i) the entanglement dynamics is 
determined only by two singlet eigenvectors of H; (ii) that 
maximal entanglement is achieved only when J rn = J[ + J' 2 ', 
(iii) the entanglement dynamics is oscillatory with period 2t*, 
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FIG. 1 : (Color online) (a) The composite system made of two sepa- 
rate Kondo spin chains initialized in their ground states in the Kondo 
regime. The extension of the clouds is tuned by J' R and J' L so that 
£fc = Nk — 1 (k — R, L). (b) To induce dynamics, one switches on 
the interaction between the two chains by the amount J m . 
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TABLE I: Typical values of J' k to generate Kondo clouds of size £ k 
as given in Eq. l[5j. 



which is only a function of J m and, thus, does not depend on 
J{ and J 2 separately. 

Generalization to a many-body system:- We now show that 
the above simple dynamics and the resulting high entangle- 
ment between the boundary spins, may be reproduced even 
with many-body systems- for arbitrary length scales- us- 
ing pertinent spin chains. For this purpose we consider two 
Kondo spin chains lllOll in the Kondo regime, i.e. two chains 
of lengths Nk described by 
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where, J\ and J 2 are nearest and next to nearest neighbor cou- 
plings, k = R (k = L) labels the right (left) chain, of is the 
vector of three Pauli operators at site i for the chain k and J' R 
( J' L ) is the impurity coupling of the right (left) hand side. 

It is well known that a Kondo spin chain supports a 
crossover from a gapless Kondo regime for J 2 < J 2 = 
0.2412 J\ to a gapped dimerized regime for J 2 > J 2 - In the 
Kondo regime the Kondo screening length is uniquely deter- 
mined by the impurity coupling Il0lll2ll and, for large chains, 
the explicit dependence is given by £j. = e"' v Jfc , where a is 
a constant; ^ sets the size of a block of spins forming a sin- 
glet with the impurity 111 211 . In the following we shall fix the 



value of J' R and J' L so that 



& = N k 



R,L. 
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We report in table I the values of the impurity couplings- de- 
termined for chains of arbitrary lengths in Ref. IU2I1 - as Nk 
is increased. Eq. (O allows to build two macroscopic singlets 
(i.e., extended over a distance £j. tuned by J' k ). The composite 
spin system is depicted in Fig. [TJ a ); the two impurities sit at 
the opposite sides and may be regarded as the boundary spins 
of the composite system while, due to Eq. (0, the two Kondo 
clouds are tuned to take over each chain separately. Note that 
not only is this J' k ~ l/Log 2 Nk much stronger than the weak 
couplings in Refs. Jj,|5|], but also the chain is gapless, so it 
cannot lead to perturbative end-to-end effective Hamiltonians. 
Initially, the two chains are separated and initialized in their 
ground states (see Fig. 01 a )) an d the initial state of the com- 
posite chain is given by IV'(O)) = Hk=R.L\GSk) where \GSk) 
is the ground state of the chain k. Then, we switch on 
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between the two chains (see Fig. |TJb)). The Hamiltonian of 
the composite system of length TV = Nl + Nr is given by 
H = Hl + Hr + Hi. Now the ground state evolves according 
to \ip{t)) = e~ lHt \ip(0)}. From knowing \ip(t)) one obtains 
the reduced density matrix of the boundary spins at a generic 
time t by tracing out all other spins from the state \ip(t)) 
and evaluate the concurrence E(t, J m ) between the bound- 
ary spins. The dynamics is now not analytically solvable and 
one has to resort to numerical simulations which, for N > 20, 
use the time dependent density matrix renormalization group 
(tDMRG) introduced in iQ while, for N < 20, one may use 
exact diagonalization. For temperatures T < l/£ ~ 2/N, 
i.e. when the two constituent chains are in the Kondo ground 
state, we find that the evolution of the composite chain well 
reproduces all the relevant features exhibited by the simple 
dynamics of a four spin system made out of two singlets. 

If the composite system built out of two extended Kondo 
singlets should reproduce the remarkable features of the sim- 
ple example discussed above one should expect that, E(t, J m ) 
oscillates with a period depending only on J m and that max- 
imal entanglement between the boundary spins is reached at 
the half of the period provided that 



J m = $(N)(J' L + J' R ). 



(7) 



$(7Y) accounts for the effects arising due to the extended size 
of the Kondo singlets. Of course, for our dynamics to make 
sense at all one has to require that, as N — s- oo, J m should 
take a nonzero and finite limiting value J (otherwise, one 
injects either zero or infinite energy). It is remarkable that 
this condition alone suffices to determine $(iV). Indeed, if 
Nl w Nr —t oo, one has 



$(ao ~ 4>g 2 (^), 
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since, in the Kondo regime, one has that £& = e" 

For the time being we will consider only chains for which 
N L = N R and so that N = 2N L - In Fig. ffla) we plot 





FIG. 2: (Color online) (a) The oscillatory dynamics of entanglement 
E vs. time t in the Kondo regime (J2 = 0) for a composite system 
of TV — 32 when Nl — Nr. (b) The maximal entanglement E m vs. 
N at t — t* when J m has its optimal value, (c) t* vs. TV. 



the evolution of the entanglement as a function of time for 
J m = 0.97Ji when N — 32 in the Kondo regime (J2 = 0) 
of each chain. We see that entanglement dynamics is oscilla- 
tory with a period 2t*. Restricting only to the first period of 
oscillations one sees that, there is an optimal value of J m for 
which, at time t*, the entanglement reaches its maximum E m , 
In Fig. |2b) we plot E m as a function of TV. Though the en- 
tanglement decreases as TV increases for short chains, its value 
remains very high and becomes distance independent for very 
long chains. It is remarkable that this distance independent 
value seems to be 0.9 (e.g. for chains of length TV = 40) 
whereas for the only other distance independent case so far 
[8[] it was merely 0.7. To complete the picture of entangle- 
ment evolution in Fig. |2jc) we plot t* as a function of TV. We 
see that the time needed to generate the entanglement between 
the boundary spins increases linearly with TV with a slope that 
is small enough to allow for fast dynamics. The linear depen- 
dence of t* on TV implies that, for a system composed of two 
extended Kondo singlets, t* is related to J m by 



FIG. 3: (Color online) (a) Optimal J m vs. TV for N L = N R in the 
Kondo regime, (b) The asymptotic behavior of <3?(TV) vs. log (-r). 



t* ~ TV ~ & 
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In Fig. [3ja) we have plotted the optimal value of J m as 
a function of TV. One sees that, as TV increases, J m goes 
to 1 thus, confirming the assumption used in the derivation 
of $(JV) (see Eq. ©). In Fig. Eft) we plot $(7V) versus 
log 2 (TV/2). The linearity of the plot provides an independent 
numerical confirmation of the result obtained in Eq. (|8). 

Our numerical approach allows to investigate also situa- 
tions for which TVj, is different from Nr. In Fig. |4{a) 
we plot E m versus TV^/TV for a composite chain of length 
N = 32. Fig. HJa) shows that the entanglement is maximal 
when Nl = Nr, and that decreases sensibly when the sizes of 
the constituent Kondo chains are very different. In Fig. Ub) 
we plot t* as a function of Nl/N. Again one sees that the op- 
timal time t* is much shorter when Nl ~ Nr. Figs. Ufa) and 




FIG. 4: (Color online) Non-symmetric case where N — 32 and Nl 
and Nr are varied: (a) E m vs. N L /N. (b) t* vs. N L /N. 



(b) lead us to conclude that efficient routing of entanglement 
is possible only if Nl ~ Nr. 

The proposed mechanism for generating high entanglement 
between the boundary spins of a composite spin system re- 
lies heavily on Eq. (O and, thus, on the fact that, for Kondo 
chains of arbitrary sizes Nl and Nr, one can always tune 
the impurity couplings J' L and J' R so as to make the Kondo 
cloud comparable with the size of the chains. As a result, 
entanglement generation between the boundary spins should 
vanish for £& < Nk/2 as well as when the constituent Kondo 
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TABLE II: Comparison between E m and t* for a Kondo spin chain 
in the Kondo ( J2 = 0) and dimer regimes (J2 = 0.42). In the table 
K stands for Kondo and D for dimer. 




FIG. 5: (Color online) A four-node router in which each user controls 
one boundary spin. A dispatcher connects two chains to induce dy- 
namics in a channel composed of two spin chains in order to generate 
entanglement between the boundary spins. 



chains are in the dimer regime (i.e., Ji > J|) where the cloud 
does not exist at all. We computed numerically E m and t*, 
for a chain composed of two Kondo spin chains in the dimer 
regime. The results are reported in table II and compared with 
the results obtained for the same quantities when the two con- 
stituent chains are in the Kondo regime and Eq. (0 is satisfied. 
Table II shows that, as N increases, entanglement E m (opti- 
mal time t*) is very small (large): for instance, for N = 40, 
in the dimer regime, E m = 0.16 and t* = 35.01 while, in the 
Kondo regime, E m = 0.89 and t* = 9.80. 

Entanglement router:- Our analysis allows to engineer an 
efficient entanglement router dispatching entanglement be- 
tween very distant qubits. A four-node entanglement router 
is sketched in Fig. [5] Each node, say A, B, C and D, has a 
boundary spin whose coupling to its adjacent chain is tuned 
so as to generate a Kondo cloud reaching the dispatch center 
(Fig. 0. The dispatcher can entangle the spins of two arbi- 
trarily chosen nodes, say A and B, by switching on a coupling 
J m between the chains A and B and, thus, induce the quench 
dynamics previously analyzed. At t = t* , the entanglement 
may be taken out of the boundary spins by a fast swap to any 
memory qubits in nodes A and B for building resources for 
quantum computation. Note that exclusive pairs of nodes, e.g. 
(A, B) and (C, D), can be connected simultaneously. 

Implementations and challenges:- Spin chains with switch- 
able/tunable couplings are realizable 11511 with both supercon- 
ducting qubits and spins in quantum dots. In the former, the 
effect of a reasonable dephasing of strength 0.005 Ji 11611 for 
N = 12 is about 10%. In the latter, a magnetic field in a 
random direction acts on each spin due to the dot nuclei 11711 . 
Our simulations show that for N = I2,a very strong magnetic 
field (~ 0.05Ji) 11711 decreases the entanglement by 5%. 

Conclusions:- We proposed a mechanism for generating 
high entanglement between distant spins by switching on an 



appropriate interaction between two Kondo spin chains. In 
contrast to other recent networking schemes [90 it does not 
demand control of the intermediate spins or time-varying lo- 
cal fields. Our results hint that a Kondo spin chain satisfying 
Eq. <(5j may be effectively described by an extended singlet 
formed by two spins since the key features of the entangle- 
ment dynamics can be easily understood using a simple model 
of a pair of two spin singlets. Indeed, in this non-solvable 
model, the above is the best explanation of the strikingly high 
entanglement. From table II one sees that in the absence of the 
Kondo cloud, entanglement is suppressed; thus, the remark- 
able dynamical behavior of the system is a new clear footprint 
of the emergence of the Kondo cloud in a Kondo system. 

Acknowledgements:- Discussions with I. Affleck, H. Johan- 
neson, N. Laflorencie and E. Sorensen are acknowledged. AB 
and SB (P.S.) thank(s) the University of Perugia (UCL) for 
hospitality and partial support. SB and AB acknowledge the 
EPSRC, and the Royal Society and the Wolfson Foundation. 



[1] S. C. Benjamin, B. W. Lovett and J. M. Smith, Laser Photon. 

Rev. 3, 556 (2009). 
[2] L. Amico, et. al., Rev. Mod. Phys. 80, 517 (2008). 
[3] A. Osterloh, et. al, Nature 416, 608-610 (2002); T. J. Osborne 

and M. A. Nielsen, Phys. Rev. A 66, 0321 10 (2002). 
[4] L. Campos Venuti, C. Degli Esposti Boschi and M. Roncaglia, 
Phys. Rev. Lett. 96, 247206 (2006); L. Campos Venuti et. al, 
Phys. Rev. A 76, 052328 (2007). 
[5] L. Campos Venuti, C. Degli Esposti Boschi and M. Roncaglia, 
Phys. Rev. Lett. 99, 060401 (2007); M. J. Hartmann, M. E. 
Reuter, M. B. Plenio, New J. Phys. 8, 94 (2006); Y. Li et. al, 
Phys. Rev. A 71, 022301 (2005); A. Wojcik et. al, Phys. Rev. 
A 72, 034303 (2005). 
[6] H. Wichterich, S. Bose, Phys. Rev. A 79, 060302(R) (2009); K. 

Sengupta, D. Sen, Phys. Rev. A 80, 032304 (2009). 
[7] F. Galve et. al, Phys. Rev. A 79, 032332 (2009). 
[8] P. Sodano, A. Bayat, S. Bose, Phys. Rev. B 81, 100412(R) 

(2010). 
[9] D. Zueco, et.al, Phys. Re v. A 80, 042303 (2009); P. J. 
Pemberton-Ross and A. Kay. I arXiv: 1 007.2786 1 
[10] E. S. Sorensen, et. al, J. Stat. Mech., P08003 (2007); N. Laflo- 
rencie, E. S. Sorensen and I. Affleck, J. Stat. Mech., P02007, 

(2008); 

[11] I. Affleck, Lecture Notes, Les Houches 2008, larXiv:0809.3474| 
[12] A. Bayat, P. Sodano and S. Bose, Phys. Rev. B 81, 064429 

(2010). 
[13] W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 
[14] S. R. White and A. E. Feiguin, Phys. Rev. Lett. 93, 076401 

(2004). 
[15] D. Giuliano and P. Sodano, Nucl. Phys. B 811, 395(2009); F. 
Meier, J. Levy, and D. Loss, Phys. Rev. Lett. 90, 047901 (2003). 
[16] J. J. G. Ripoll, et. al, Phys. Rev. B 77, 024522 (2008). 
[17] J. M. Taylor, et. al, Phys. Rev. B 76, 035315 (2007). 



